arXiv:l506.04164v2 [gr-qc] 9Jul2015 


Noncommutative spaces and covariant formulation of statistical mechanics 

V. HosseinzadehJB M. A. GorjiQ and K. NozarH 
Department of Physics, Faculty of Basic Sciences, 

University of Mazandaran, P.O. Box f7fl6-954f7 Babolsar, Iran 

B. Vakil§ 

Department of Physics, Central Tehran Branch, Islamic Azad University, Tehran, Iran 

We study the statistical mechanics of a general Hamiltonian system in the context of symplectic 
structure of the corresponding phase space. This covariant formalism reveals some interesting cor¬ 
respondences between properties of the phase space and the associated statistical physics. While 
topology, as a global property, turns out to be related to the total number of microstates, the in¬ 
variant measure which assigns a priori probability distribution over the microstates, is determined 
by the local form of the symplectic structure. As an example of a model for which the phase space 
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I. INTRODUCTION 

Statistical mechanics is a bridge between mechanics 
and thermodynamics. Indeed, it may be assumed as the 
microscopic mechanical basis for the macroscopic ther¬ 
modynamical properties of a physical system. The key 
concept which links these two arenas is the notion of 
the phase space, that is the space of all distinct pos¬ 
sible microstates. The dynamical evolution is nothing 
but a smooth transition from one microstate to another 
and statistical mechanics is about how to count the num¬ 
ber of these microstates. In a more technical language, 
the phase space of a dynamical system is a symplectic 
manifold which naturally is equipped with a symplectic 
structure .1]. The symplectic structure covariantly deter¬ 
mines the Poisson brackets and Liouville volume and the 
kinematics is then defined in a coordinate-independent 
manner on the phase space. The dynamics will be also 
specified upon taking the Hamiltonian function as the 
generator of time evolution in this setup. Therefore, the 
Hamiltonian determines how the system evolves through 
the microstates and the Liouville volume specifies how 
one can count and measure them. Consider, for in¬ 
stance, a mechanical system consisting of N particles 
which are subjected to some forces. Each particle moves 
in Euclidean three-dimensional space R 3 and the corre¬ 
sponding configuration space is then R 3Ar . Therefore, the 
phase space will be the cotangent bundle of this space, 
that is, H 3N x R 3Ar , which naturally admits a symplec¬ 
tic structure. In terms of the positions and momenta of 
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the particles, the associated symplectic structure takes 
the standard canonical form. Such a canonical symplec¬ 
tic structure, which is also dynamically invariant, may 
be used to determine the fundamental canonical Poisson 
brackets and the kinematics for the system under consid¬ 
eration. Subsequently, based on the canonical symplectic 
structure one leads to the canonical measure which can 
be viewed as a mathematical expression for the principle 
of equal probabilities in standard statistical mechanics 
0. Although it might be thought that the symplectic 
formulation is nothing but a covariant reformulation of 
the standard statistical mechanics, its advantages may 
become more clear when one is dealing with the Hamil¬ 
tonian systems with nontrivial structures (for instance 
when the standard R" topology of the phase space is re¬ 
placed with a more complicated topology or the standard 
canonical Poisson algebra is replaced with a deformed 
noncanonical one). Such systems, are investigated in the 
context of doubly special relativity theories, as a flat limit 
of quantum gravity, where the four-momentum space has 
a curved (de Sitter or anti-de Sitter) geometry with non¬ 
trivial topology [3]. The space-time structure then turns 
out to be noncommutative and the associated Poisson al¬ 
gebra between the positions and momenta become non¬ 
canonical SI- The space-time with noncommutative 
coordinates was first proposed by Snyder in Ref. Q in 
order to regularize quantum field theory in the ultravio¬ 
let regime pj. In more recent times, it was also used by 
the stability theory of the Lie algebras to be a minimal 
candidate for the relativistic quantum mechanics Q and 
by string theory at a fundamental level |§]. The duality 
between curved momentum space and noncommutative 
structure was first pointed out by Majid in Ref. 0- The 
direct consequence of space-time with noncommutative 
coordinates is the modification to the dispersion relation 
and the associated density of states [lH - fla ] . In symplec- 
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tic formulation, the density of states is properly deter¬ 
mined by the Liouville measure which will be constructed 
from the symplectic structure H3- In noncommutative 
space-time, however, the corresponding symplectic struc¬ 
ture seems to be noncanonical which induces a nonuni¬ 
form probability distribution over the set of microstates 
through a deformed Liouville volume element. This dy¬ 
namically invariant Liouville volume enters in the defi¬ 
nition of the Gibbs entropy and partition function from 
which all the thermodynamical properties of a system can 
be extracted. Therefore, the topology of the phase space 
will play an important roll since all integrals are eval¬ 
uated over the entire phase space. For instance, let us 
consider a system which has an ultraviolet cutoff due to 
the existence of a maximal momentum. It is easy to see 
that the topology of the spatial sector of the momentum 
space may be compact. Moreover, systems with a totally 
compact phase space, which can be considered as the 
classical limit of quantum systems with finite dimensional 
Hilbert space (e.g, angular momentum for systems with 
fixed total angular momentum) were recently studied in 
the context of loop quantum gravity [18] . For such a com¬ 
pact phase space, the Liouville volume is finite and con¬ 
sequently a finite number of microstates naturally arises 
jl9l ]. This shows the correspondence between the topol¬ 
ogy of the phase space and the number of microstates 
which can be realized when the symplectic geometry is 
implemented. The black hole physics is a good exam¬ 
ple in which the system obeys the thermodynamical laws 
and its statistical origin is not thoroughly formulated yet 
sain in the framework of Hamiltonian formalism (see 
also Ref. @)- 


In this regard, formulating the statistical mechanics in 
its most fundamental form is important, at least, for two 
reasons: (i) finding a more precise and fundamental inter¬ 
pretation for the basis of the statistical mechanics, and 
(ii) studying the statistical mechanics for the Hamilto¬ 
nian systems with nontrivial structure, such as those we 
have mentioned above. Motivated by the stated issues, 
in this paper we are going to formulate statistical me¬ 
chanics for a general Hamiltonian system in a covariant 
(coordinate-independent) manner. The paper is orga¬ 
nized as follows. In section II, we deal with the kinemat¬ 
ics and dynamics of a many-particle system with general 
Hamiltonian structure in the context of symplectic geom¬ 
etry. The statistical mechanics for such systems is formu¬ 
lated in section III. In section IV, by means of the con¬ 
structed setup we study the statistical mechanics in the 
Snyder noncommutative space-time with curved energy- 
momentum space as an explicit example of a physical 
system with nontrivial (noncommutative) Hamiltonian 
structure. As a case study, the thermodynamical prop¬ 
erties of a harmonic oscillator system in Snyder space is 
also presented at the end of this section. Section V is 
devoted to the summary and conclusions. 


II. HAMILTONIAN SYSTEMS 

In this section, in order to study the statistical mechan¬ 
ics in the context of symplectic geometry, we consider the 
kinematics and dynamics of an V-particle system. 

A. Kinematics 

Let us consider a system of N particles similar to con¬ 
ventional systems in statistical mechanics. The corre¬ 
sponding 6./V-dimensional kinematical phase space may 
be obtained by coupling N single-particle phase spaces 
r a with a = 1,..., N as 

T = T 1 x...xT n . (1) 

The phase space T naturally admits a symplectic struc¬ 
ture uj which is a nondegenerate closed 2-form [lj. The 
classical state of this V-particle system is determined by 
a point in T (phase point). The evolution of phase points 
is then determined by the Hamiltonian vector field x H , 
whose integral curves are trajectories of the phase points 
(phase trajectories). From the fact that w is nondegen¬ 
erate, one can assign a vector field to a function / as 
w(x f ) = df and the Poisson bracket between two observ¬ 
ables (real-valued functions on the phase space) can be 
defined as 

{/, 9} =w(x / ,xj. (2) 

These functions with the above structure constitute a Lie 
algebra under Poisson brackets on T. We also have the so- 
called Liouville volume on T with its standard definition 
as 

uj 3N = uj A ... A uj (3 N times). (3) 

In terms of positions and momenta x = (q, p) of the 
particles, the Liouville volume takes the local form 

uj 3N = Vdet uj dq\ A ... A dq 3 N A dp\ A ... A dp% = 

Vdetu >d 3N qd 3N p, (4) 

where q l a and pf are the Pth component of the positions 
and momenta of the crth particle, respectively, and u> is 
the matrix representation of the symplectic structure uj 
with components u,j = uj(d/dx l a ,d/dx J p). 

Now, one may be able to decompose the symplectic 
structure as 

N 

w = (5) 

a=l 

where uj a is the symplectic structure on r a . Although 
this assumption has no mathematical basis, from a phys¬ 
ical point of view it is acceptable since the particles are 
assumed to be kinematically separated. In other words, 
the Lie algebras of the functions defined by {uj a } over 
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{T„} are separately closed. Therefore, the Liouville vol¬ 
ume over T takes the form 

<6) 

where the wedge products take place 3 N times. The 
volume element then works out to be 

co 3N = ujf A ... Acoff = T JJ \/det d 3N q d 3N p , (7) 

where w 3 is the volume element of the corresponding 
{T q }, and we have also used the fact that w™ = 0, for 
n > 3. Clearly, Vdetcu, that is, a function on T, factor¬ 
izes into products of {\/det ui a }. This cannot be realized 
in the most general situation, where the symplectic struc¬ 
ture cannot be written as ®. 

As a special case, consider a system consisting of N 
particles subject to some forces moving in Euclidean 
three-dimensional space with standard R 3 topology. 
Therefore, the phase space will be R 3jY x R :i :V which is 
indeed the cotangent bundle of the configuration space, 
that is, R ;V . Geometrically, the cotangent bundle is a 
symplectic manifold endowed with a canonical symplectic 
structure 


w c 


dQ' a A dP° 


( 8 ) 


in which Q and P are interpreted as the positions and mo¬ 
menta of particles, respectively. The summation is over 
both i = 1,2,3 and a = 1, This is the canonical 

representation of symplectic structure and coordinates 
(QL,Pf) are known as the canonical coordinates. Sub¬ 
stituting the canonical structure © into the definition 
©, the Poisson bracket will be 


df dg df dg 

dQL <9P“ <9P ? dQi ’ 


( 9 ) 


which leads to the standard canonical Poisson algebra 


{Qa.QJ}c=o, {Q^h = {Pf,P 3 p } c = o. 


( 10 ) 


Also, with the help of © we get from © the Liouville 
volume as 

ui 3N = d 3N Q d 3jV P , (11) 

which is nothing but the standard volume element on 
R 6JV . Thus, the kinematics of the standard statistical 
mechanics can be recovered as a special case of this setup. 
Taking into account the fact that det u) c = 1, the measure 
associated to the standard volume element m assigns a 
uniform probability distribution to the set of microstates. 
This is justifiable in the light of Laplace’s principle of in¬ 
difference, which states that in the absence of any further 
information, all outcomes are equally likely pj. This is 
the fundamental basis of the statistical mechanics. 


In the case where the space-time has a noncommu- 
tative structure, apart from the details of the different 
models, there is always a deformation parameter which 
can be linked to a minimal length associated with the 
system under consideration. The direct consequence of 
this setup is the modification of the dispersion relation 
and the density of states (see Ref. [8j for the special case 
of the stable noncommutative algebra for the relativis¬ 
tic statistical mechanics). The density of states is de¬ 
termined by the symplectic structure and the associated 
Liouville volume. Since the symplectic structure takes 
a noncanonical form in terms of the physical positions 
and momenta (q,p) in noncommutative spaces [l9f, we 
have det a; = det uj(q,p) ^ 1. This result immediately 
shows that Laplace’s indifference principle is no longer 
applicable for a noncommutative space-time. This is be¬ 
cause of the fact that there is extra information in these 
setups (minimal length or maximal momentum) and the 
Liouville measure then assigns a nonuniform probabil¬ 
ity distribution over the set of microstates by means of 
the relation ©■ This simple consideration of noncom¬ 
mutative space-time as an example of a kinematically 
deformed system shows the advantage of the coordinate- 
independent symplectic formulation of the statistical me¬ 
chanics (see Refs. [UHII for explicit examples). We will 
explicitly consider such an example in section IV. 

It also should be noted that, according to the Dar- 
boux theorem, there is always a local chart on T in which 
any symplectic structure takes the canonical form with 
detu> = 1. In the case of kinematically deformed sys¬ 
tems such as the noncommutative space-time, however, 
the new canonical coordinates cannot be interpreted as 
the positions and momenta of the particles. Moreover, in 
the case of standard statistical mechanics, one can also 
find a local chart in which the standard canonical struc¬ 
ture w c takes a noncanonical form with detu; c ^ 1. But 
again, these new noncanonical coordinates cannot be in¬ 
terpreted as positions and momenta of particles. In all 
of the above cases, as we will see in the next section, the 
partition function and consequently resultant thermody¬ 
namical properties are independent of the chart in which 
the system is considered. 


B. Dynamics 

As we have mentioned above, the time evolution of the 
system is determined by the Hamiltonian vector field x H . 
It can be obtained by solving the dynamical equation 

oj (x H ) = dH , (12) 

where H is the Hamiltonian function of the system that 
is a real-valued function over T and w (x H ) is an interior 
product of uj and x H . Additionally, x H can be expanded 
as 

N 

x h = 5I x S’ ( l3 ) 

a—1 
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where x“ is the projection of x H on T a . Note that, in 
general, the domain of components of x“ is the whole 
phase space T and consequently the Lie brackets between 
{x“} will not vanish: 

o. (14) 

However, as a special case, one can consider a system con¬ 
sisting of noninteracting particles. Therefore, the Hamil¬ 
tonian function can be written as the sum of individual 
Hamiltonians, that is, 

N 

H = Y,H a - (15) 

c*=l 

So, from decomposition © for the symplectic structure, 
the dynamical equation m becomes 

(L w “)(E x ») =E dHo ■ (16) 

'c*=i / k /3=1 / a=1 

Using the fact that w Q (x^) = 0, for a/|3, one is led to 
a set of N independent differential equations 

wi(x^) = dffi, ... , wjv(x^) = dH N . (17) 

Now, it is clear that the domain of the components of 
x“ will be T a and the Lie brackets between x“ and x() 
vanish 

K>x£] = 0. (18) 

This can be seen as a criterion for the kinematical and 
dynamical separability of the particles. In the Appendix, 
we will see that this benchmark shows itself as a crite¬ 
rion for the statistical independence of the macroscopic 
subsystems which in turn is a definition of equilibrium. 

The other important property of the Hamiltonian sys¬ 
tems is the existence of an invariant measure that en¬ 
ables one to construct the equilibrium statistical mechan¬ 
ics through the well-known Liouville theorem as 

5 “” = °’ < 19 > 

where d/dt = d/dt + C x is the total time derivative and 
£ Xff denotes the Lie derivative with respect to x H . The 
conservation of ui 3N can be traced back to the fact that 
C Xh oj 3N = 3 N(C Xh w) A lu 3N ~ 1 and = du (x H ) + 

d(u (x H )) = 0, where we have used d(w (x a )) = di 2 H = 0 
and the closure of the symplectic structure dco = 0. It is 
obvious that oj 3N is not explicitly time dependent. Note 
also that the Liouville theorem arises from the particular 
dynamics of a Hamiltonian system. 

Consider the standard statistical mechanics as a spe¬ 
cial case, where the kinematics is given by canonical sym- 
plectic structure ©. The dynamical equation (fT2l) then 
leads to the familiar form of the Hamilton’s equations in 
terms of the physical positions and momenta as 

d(& = dH_ dPl = _dH_ 

dt dPf ’ dt dQ\ 1 ’ 


Note that when the symplectic structure has a noncanon- 
ical form, the form of Hamilton’s equations would deviate 
from the relation (I20|) . 

The triple (T,oj,H) constitutes a general Hamiltonian 
system and we explore the statistical physics of this gen¬ 
eral system in the next section by means of an invariant 
volume element oj 3N which is covariantly defined in the 
relation ©■ 

III. STATISTICAL MECHANICS 

Statistical mechanics links macroscopic properties of a 
system to microscopic laws that are quantum mechanical 
at the fundamental level. The basic concept is the von 
Neumann entropy that is defined as 

S=-Tr(p\np), (21) 

where p is the density operator [24j |. Using the princi¬ 
ple of maximum entropy and considering specific relevant 
statistical constraints, we can find the density operator 
which contains all the statistical information about the 
system at equilibrium. For example, for a system in con¬ 
tact with a thermal bath at temperature T and therefore 
with a fixed mean energy, the resultant normalized den¬ 
sity operator works out to be 

P* = \ exp [-H/T ], (22) 

where H is the Hamiltonian operator of the system and Z 
is called the canonical partition function which is defined 
as 

Z = Tr exp [~H/T\ = ^ exp [ - Ei /T] , (23) 

i 

where {£j} are the energy eigenvalues of H and the sum¬ 
mation is taken over all the accessible microstates for 
the system. At the fundamental level, the microscopic 
laws are quantum mechanical and the microstates are 
determined by quantum mechanics. On the other hand, 
in the classical limit, the microstate of a system repre¬ 
sents itself as a point on the corresponding phase space 
and due to the infinite resolution of phase space points, 
the summation over microstates is really not well defined. 
Consequently, one cannot construct a full classical statis¬ 
tical mechanics in essence. However, quantum mechanics 
provides the semiclassical approximation in the spirit of 
the uncertainty principle. Indeed, a phase space with 
finite resolution and well-defined number of classical mi- 
crostates can be achieved through the approximation E3 

<24) 

where 1 /(AT!) is due to the fact that the particles may 
be indistinguishable and ui 3N is the dynamically invari¬ 
ant Liouville volume element of the 61V-dimensional sym¬ 
plectic manifold T which is defined in the relation © . It 
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is important to note that the semiclassical approxima¬ 
tion (1241) coincides with the full quantum consideration 
at high temperature regime. In this limit, the von Neu¬ 
mann entropy m is replaced with the Gibbs entropy 

s = (25) 

where p is now interpreted as a probability density over 
the space of microstates T. The maximum entropy prin¬ 
ciple for a system in heat bath at temperature T leads to 
the Gibbs state 

Pc = Tf exp [~H/T \, (26) 

where H is the Hamiltonian function of the system and 
Z is the semiclassical canonical partition function 

Z = M 6XP [ ~ H / T ]u 3N • (27) 

The result of the integral over the phase space for the 
total partition function m is independent of a chart 
in which the system is considered when the dynamical 
equation m is satisfied in a chart-independent man¬ 
ner. While the local form of the symplectic structure 
determines the probability distribution over the set of 
microstates, from the statistical point of view, we expect 
that the resultant thermodynamical properties extracted 
from a partition function are independent of a chart in 
which the physical system is considered (see, for instance, 
Ref. G3 in which the partition function in two different 
charts is calculated). 

By considering the decompositions ([5]) and (fl5l) for the 
symplectic structure ui and Hamiltonian function H , re¬ 
spectively, the probability density will also be decom¬ 
posed as 

p c ^ = (p^?)A...A(pf^), (28) 

where is the Gibbs state defined on T a . Hence the 
integral over T can be factorized into the products of the 
integrals over the phase spaces of the particles as 

z m n (l ex p [ - h ^ t ] <4) • (29) 

Now, if the particles are influenced by the same kinemat¬ 
ics and dynamics, the partition function can be written 
in the well-known form 


where we have defined the single-particle state partition 
function 

Zi= [ eM-Hi/T] w? . (31) 

Jr 1 

It is usually claimed that the relation (13U1) is applicable 
when one ignores the quantum correlations between par¬ 
ticles. This means that decompositions such as © and 


(fl5l) may be considered as conditions for the absence of 
quantum correlations. Moreover, the general formalism 
which is formulated in this section, can also support a 
Hamiltonian system with very nontrivial structure, e.g. 
a kinematically classical entangled system in which the 
decomposition © is no longer applicable. 


IV. APPLICATION IN NONCOMMUTATIVE 
SPACES 

Existence of a minimal length is suggested by any 
quantum gravity candidate such as string theory and loop 
quantum gravity [2(| . It is then widely believed that 

a nongravitational theory which supports the existence 
of a minimal length would arise at the flat limit of quan¬ 
tum gravity. Such a theory, would be reduced to the 
standard relativistic quantum mechanics at low energy 
regime where the effects of a minimal length are negli¬ 
gible. Evidently, taking a minimal length into account 
naturally leads to the space-time with noncommutative 
coordinates. The first attempt in this direction was done 
by Snyder in 1947 who formulated a Lorentz-invariant 
discrete space-time with noncommutative coordinates @. 
The space-time with noncommutative structure is also 
suggested by the theory of stability of the Lie algebras 
0. On the other hand, existence of a universal mini¬ 
mal length cannot be supported by the special relativity 
since any length scale in one inertial frame may be differ¬ 
ent in another observer’s frame through the well-known 
Lorentz-Fitzgerald contraction. Thus, the doubly spe¬ 
cial relativity theories are formulated in order to take 
into account an observer-independent minimum length 
as well as the velocity of light [27| . The curved four- 
momentum space then naturally appears to be a suitable 
framework to formulate the doubly special relativity the¬ 
ories [3]. Furthermore, the different doubly special rel¬ 
ativity theories can be obtained from the different basis 
of the K.-Poincare algebra on the associated K-Minkowski 
noncommutative space-time [28| . In this respect, it is also 
shown that the Snyder noncommutative algebra which 
is proposed in Ref. @ and the stable algebra of rela¬ 
tivistic quantum mechanics that is obtained in Ref. jj 
can be obtained from a ten-dimensional phase space with 
curved geometry for the four-momentum space through 
the symplectic reduction process Q. Thus, it seems that 
the space-time with curved four-momentum space and 
noncommutative coordinates is a fundamental framework 
for taking into account a minimal length scale in the flat 
limit of quantum gravity [29}. It would also be claimed 
that the Lorentz invariance is an approximate symmetry 
and will be broken at high energy regime. Therefore, the 
Lorentz-violating noncommutative algebra was proposed 
by Camelia in Ref. [27j which supports the existence of 
a minimal observer-independent length scale. Magueijo 
and Smolin showed that the Lorentz invariance can also 
be preserved by a nonlinear action of the Lorentz group 
on the momentum space (30). In order to preserve the 
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Lorentz invariance, the following commutation relations 
for the generators will be held 

{Jab: Jcd\ — JadPbc T JbcVad JbdPac Jac^bd: (32) 

{Jab: Pc } — PaPbc Pbdac: {Jab: X c } — Xa^Jbc XbVac: 

where r] a b = diag(+l, —1, —1, —1) is the Minkowski met¬ 
ric and a, &, c, d = 0, The commutation rela¬ 

tions between ( x a ,p a ) then classify the different Lorentz- 
invariant_algebras such as the Snyder noncommutative 
algebra fb], the stable algebra [3| and also the standard 
Lorentz algebra. 

While the relativistic algebras are written in the eight¬ 
dimensional relativistic phase space, to formulate the sta¬ 
tistical mechanics one needs to work in a six-dimensional 
nonrelativistic subalgebra of the deformed relativistic al¬ 
gebras in which the time parameter is fixed and the cor¬ 
responding Poisson brackets take a noncanonical form. 
By means of the resulting deformed Hamiltonian system, 
one can study the statistical mechanics with the help of 
the constructed setup in the pervious sections. In recent 
years, many works have gone in this direction in order to 
study the effects of an ultraviolet cutoff (minimal length 
and maximal momentum) on the thermodynamical prop¬ 
erties of the physical systems. For instance, thermo¬ 
dynamics of some physical systems in noncommutative 
spaces is studied in Ref. & For the special case of the 
Snyder noncommutative space, see Ref. [32| . in which 
the thermodynamics of the early Universe is explored. 
It is shown that the Liouville theorem for the deformed 
phase space with general noncanonical Poisson algebra is 
satisfied. However, as we have shown in section III, the 
Liouville theorem can be justified in a very simple manner 
in the covariant formalism, which also shows the advan¬ 
tage of the setup. Furthermore, thermodynamical ideal 
gases are studied in the context of doubly special rela¬ 
tivity theories [16]. Also, motivated by the string theory, 
the generalized uncertainty principle is suggested which 
supports the existence of a minimal length as a nonzero 
uncertainty in position measurement |25l |33j[ . For the 
statistical mechanics in this setup, see Refs, fuilrt. In¬ 
spired by loop quantum gravity, the polymer quantum 
mechanics is formulated which supports the existence_of 
a minimal length known as the polymer length scale [36( |. 
Thermodynamical properties of the ideal gases and har¬ 
monic oscillator are also studied in Refs, [F3, |37| . To see 
the effects of minimal length on the thermodynamics of 
the black holes in noncommutative space, the generalized 
uncertainty principle framework, and the polymer quan¬ 
tization scheme, see Refs. [38l-(40|. Apart from the details 
of the above-mentioned effective approaches to quantum 
gravity phenomenology, all of them try to consider the 
effects of a minimal length in a relevant manner, and im¬ 
plementing the covariant formalism can clarify the con¬ 
sequences and applications in its most fundamental form 
in the language of symplectic geometry. 


A. Statistical mechanics in the Snyder 
noncommutative space 

As we have mentioned above, the advantages of the co¬ 
variant formulation may be revealed when a Hamiltonian 
system with nontrivial structure is considered. There¬ 
fore, in this subsection we study the statistical mechanics 
in the Snyder noncommutative space as a well-known ex¬ 
ample of a kinematically deformed Hamiltonian system. 

The Snyder relativistic algebra preserves the Lorentz 
invariance and therefore the commutation relations (1321) 
are satisfied by the corresponding generators. The com¬ 
mutation relations between (x a ,p a ) which define the Sny¬ 
der algebra are given by f(i 

{Xa:X b }= {x a ,p b } = TJab + {Pa:Pb}= 0, 

Av Av 

(33) 

where k is the deformation parameter with the dimen¬ 
sion of the inverse of length which plays the role of 
the universal quantum gravity scale in this setup. The 
four-momentum space of the eight-dimensional relativis¬ 
tic phase space of a test particle which moves on the 
space-time with noncommutative algebra (1331) . is a de 
Sitter space with topology RxS 3 [!, la] • Identifying the 
energy space with R, the topology for the space of spa¬ 
tial momenta will be S 3 . In order to study the statistical 
mechanics in this setup, one should consider the non¬ 
relativistic subalgebra of (1351) and replace the standard 
canonical Poisson algebra with them. The nonrelativistic 
Snyder algebra for a particle moves in three-dimensional 
Euclidean space R 3 and is then given by (4l| 

J ■ ■ VV' 

{QiiQj} = 2 ^ ’ Pj 1 = ~2~ 1 i Pj } = 0* 

AC AC 

(34) 

It is straightforward to show that the above commutation 
relations can be realized from the symplectic structure 
0 

U) = dq l A dpi — — d{q l pi) A d In (p 2 + k 2 ) , (35) 

through the covariant definition ([5]) for the Poisson brack¬ 
ets. The associated phase space of the particle is then 
simply the cotangent bundle of the configuration space 
which now has the nontrivial R 3 x S 3 topology. In¬ 
deed, the space of the spatial momenta has compact 
S 3 topology rather than the standard R 3 one (see Refs. 
Pj3>EH,[hi] for more details). Since the quantum gravity 
cutoff (which is defined by the deformation parameter k) 
is universal, it will be the same for all of the particles. 
Thus, for the kinematically deformed system consisting 
of N particles which obeys the Snyder noncommutative 
algebra ED, the total symplectic structure can be ob¬ 
tained by the relation The phase space r a for all of 
the particles is the same and has the nontrivial R 3 x S 3 
topology. The associated total phase space Ts can be 
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easily obtained through the coupling CD, which has the 
following nontrivial topology: 

R 3N x S 3 x ... x S 3 . (36) 

N times 


Substituting the symplectic structure ( 1551 ) into the defi¬ 
nition ©, the corresponding Liouville volume takes the 
form 


>3 n = TT ( d 3 q a A d 3 p a \ = d 3N q A d 3N p 

Ji v (i + (p %) 2 )) nLi (i + (p“A) 2 ) ’ 

(37) 


where q l a and pf are the i’th component of the posi¬ 
tions and momenta of the ath particle respectively. From 
this relation it is clear that the probability distribution 
is nonuniform in Snyder noncommutative space and the 
standard uniform one can be recovered only at low en¬ 
ergy (low temperature) regime k —> oo. Indeed, the 
existence of a minimal length, as an extra information, 
changes the probability distribution at the high energy 
regime. Thus, in all of the kinematically deformed non¬ 
commutative phase spaces, the local form of the symplec¬ 
tic structure (in terms of the positions and momenta of 
the particles) is noncanonical and therefore the proba¬ 
bility distribution will be nonuniform such that the mi¬ 
crostates with higher momenta are less probable. This is 
the advantage of the symplectic covariant formulation of 
the statistical mechanics which shows that the probabil¬ 
ity distribution will be nonuniform on any phase space 
with noncanonical Poisson algebra like the Snyder alge¬ 
bra ( 1551 ) . However, in light of the Darboux theorem, one 
can always find a local chart through the Darboux trans- 


qjyPi 


= 


formation (g£,pf) -t (X z a = q^ - {pa) * +K * » y k , 
pf), in which the symplectic structure ( 1551 ) takes the 
canonical form = dX^AdYf* for all of the parti¬ 
cles. Thus, the associated Liouville volume becomes 
to 3N = d 3N X A d 3N Y, which induces a uniform prob¬ 
ability distribution over the set of microstates. But, it 
is important to note that (X l a ,YA) are different from 
the standard canonical ones (Q l a ,P!f) which are defined 
in relations © and (fill) . Indeed, while (Q),,P“) are 
the positions and momenta of the particles by definition, 
(X l a , Yf) cannot be interpreted as the positions and mo¬ 
menta of the particles. Also, the noncanonical Snyder Li¬ 
ouville volume (1371) reduces to the standard canonical one 
m in the limit of n —> oo and the noncanonical variables 
(q l a ,pf) then coincide with standard canonical variables 
(Q^P") in this limit. However, the canonical variables 
(X^,IA“) are obtained through a Darboux transforma¬ 
tion and the functional form of the Hamiltonian function 
will be modified (compared with the standard functional 
form) in this chart as H(X, Y) since the dynamical equa¬ 
tion m will be satisfied in a chart-independent manner. 

Substituting the Liouville volume ( 1571 ) in the relation 


the total partition function will be 


Z 


1 

N\ 


r exp [~H(q,p)/T] 

k nli (i + 0p%) 2 ) 


d 3N qd 3N p. 


(38) 


Since the Snyder Liouville volume ( 1571 ) is decomposed 
similar to the relation CD, for the noninteracting sys¬ 
tems in which the total Hamiltonian function also de¬ 
composes as the relation CSD, one can rewrite the to¬ 
tal partition function ( 1551 ) in terms of the single-particle 
partition function. Taking the decomposition m into 
account and substituting the volume ( 1571 ) into relation 
(1381) . the total partition function can be rewritten in the 
form of the relation ( 1551 ) as 



exp [ - H a (q a ,p a )/T] 3 3 


[l + (p“/k ) 2 


■ d q a d p c 


(39) 


where we have defined the single-particle partition func¬ 
tion as 


Z! 



f ,3 i exp [ — Hi(qi,p 1 )/T] 
J S3 P (l + (p 1 /^) 2 ) 


(40) 


and we have also assumed that the Hamiltonian func¬ 
tions are the same for all of the particles. Then, all the 
thermodynamical properties of a system can be deduced 
from the total partition function ( 1551 ) or ( 1551 ) in which, as 
we have mentioned in the pervious section, the Gibbs fac¬ 
tor 1/N\ is considered for the nonlocalized systems such 
as the ideal gas and it should then be removed for the 
localized systems such as the harmonic oscillator which 
is the subject of the next subsection. 


B. Thermodynamics of 3D harmonic oscillator 


Now let us consider a system of N independent and 

similar three-dimensional harmonic oscillators. Since the 

oscillators are assumed to be independent, the relations 

© and (1151) will be satisfied and one then can implement 

the relation (1391) in order to obtain the partition function 

of the system. The Hamiltonian of the three-dimensional 

2 

simple harmonic oscillator is given by H\(q,p) = + 

i ma 2 q 2 , where m is the mass of the oscillator, a is the 
frequency, p 2 = p 2 + p 2 + p 2 and q 2 = x 2 + y 2 + z 2 . 
Substituting this Hamiltonian into the relation (1401) , gives 
the following single-particle partition function 




2 r j ~<‘2 

K 1 

mh 3 cr 3 




yj2mT 


erfc 


K 

\/2 mT 


exp 


2mT 



At the high temperature regime, the second term in the 
right-hand side of the above relation is negligible. Under 
this condition the dominant contribution is due to the 



















first term and the partition function is approximated by 
Z 1 « 3 ■ This shows that in the domain of the high 

temperature, the number of degrees of freedom will be 
reduced since the partition function is proportional to 
~ T 2 rather than the standard one which is proportional 
to ~ T 3 . We will show this feature in a more precise 
manner in this subsection. The total partition function 
can be obtained from the relation (l39l) as 

^ = Zi N , (42) 


where the Gibbs factor 1/N\ is dropped since the os¬ 
cillators are localized. Substituting the single-particle 
partition function (1411) into this relation, the total parti¬ 
tion function for N three-dimensional harmonic oscilla¬ 
tors can be obtained from which one is able to study the 
thermodynamical behavior of the system in Snyder space. 
However, before doing this task, it is useful to consider 


the limit 


1 which is related to the first order cor¬ 


rections that would arise from the minimal length effects. 
This limit also allows us to compare our result with the 
full quantum partition function. In this limit, we may use 


the relation erfc(ar) = E^L 0 (~ 1 ) n ( ( to arrive 

at the following asymptotic expression for the partition 
function 


Zi 



71=1 


(2n+l)!l\ 

( 20 2 )™ J ’ 


(43) 


in which we have defined the dimensionless variable 0 = 
-t== . In order to obtain the quantum partition func¬ 
tion, one should solve the Hamiltonian eigenvalue prob¬ 
lem in Snyder space which is not an easy task at all since 
the Schrodinger equation takes a complicated form in this 
setup. However, in the context of the generalized uncer¬ 
tainty principle the problem of a Z?-dimensional harmonic 
is exactly solved j44| and the results are applicable also 
for the Snyder algebra flHD as a particular case. The 
energy eigenvalues for a three-dimensional harmonic os¬ 
cillator in Snyder space are given by 


E n i = ho 


17 3\ 

\A+ 

/ mho 

\ 2 

.(”+ 2 ) 

l 2k 2 

) 

mho / 
+ 2k 2 ( 

(n + 

-Y- 

2 ) 

0( /+i )+B) 


(44) 


u 

N 



Non-deformed 

Snyder-deformed 


FIG. 1. Internal energy versus temperature is plotted. The 
solid and dashed lines represent the internal energy for the 
Snyder-deformed and nondeformed cases, respectively. The 
deviation from the standard case arises at the high tempera¬ 
ture regime when the minimal length effects dominate. The 
figure is plotted for m = 1 = k. 


it leads to Z\ = (^) 3 (l — §0 -2 ) + O [0 -4 ], which 
is nothing other than the first term in the summation 
(1431) . This coincidence shows that the quantum partition 
function (H5l) correctly leads to the semiclassical parti¬ 
tion function m up to the first order of approximation 
at the high temperature regime. It also reveals the ad¬ 
vantages of semiclassical approximation (HH1) in the sense 
that although there is not an analytical solution for the 
case of full quantum partition function (1451) . the rela¬ 
tion (1411) provides an analytical expression for the high 
temperatures. Therefore, in dealing with the statistical 
considerations of the physical systems in the presence of 
a minimal length, since one is usually interested in the 
high temperature regime (such as early Universe thermo¬ 
dynamics), the usage of the semiclassical approximation 
seems to be quite reasonable. 

Substituting the single-particle partition function aud 
into the relation (TT21) gives the total partition function. 
The internal energy then can be obtained from the stan¬ 
dard definition U = T 2 N as 


U = 


3 NT 


1 + - ^1 — \ZttQ erfc[0] e e ^ — 0 


.(46) 


where n = 0,1,2,..., and 0 < l < n and also we have 
considered relevant parameter identification to rewrite 
the result of the Ref. 0] in our notation and units. 
The associated quantum partition function then will be 
obtained from the definition (E^l) as 


z i = EE ex p 

n=0 l —0 


Enl 

~T~ 


(45) 


The above summation cannot be evaluated analytically. 
However, in Ref. [35;. it is shown that up to the first or¬ 
der of approximation, in the limit of 0 = ^== 1, 


In figure [T| we have plotted the internal energy versus 
temperature in comparison with its nondeformed coun¬ 
terpart. The specific heat can also be deduced from the 
definition C v = (jpp) v which leads to a somehow cum¬ 
bersome expression. 

An interesting feature here is the reduction of the num¬ 
ber of degrees of freedom in Snyder space. According 
to the well-known equipartition theorem of energy, for 
the Hamiltonians of the form Hi(q,p) = ^ ma 2 q 2 , 

each of the 6 degrees of freedom makes a contribution 
of \T towards the internal energy. Thus one may con¬ 
sider as the number of degrees of freedom for each 

























9 


2 U 
NT 



FIG. 2. The number of degrees of freedom versus tempera¬ 
ture. As it is clear from the figure, in Snyder space the number 
of degrees of freedom for a three-dimensional harmonic oscil¬ 
lator is reduced from 6 to 4 at the high temperature regime. 


three-dimensional harmonic oscillator. In figure [2] the 
number of degrees of freedom for each three-dimensional 
harmonic oscillator is plotted. As this figure shows the 
number of degrees of freedom will be reduced from 6 to 
4 in Snyder space. This may be compared with the sta¬ 
tistical mechanics of the ideal gases considered in Ref. 
phi ] in which the same result is obtained. However, we 
would like to emphasize that the reduction of the num¬ 
ber of degrees of freedom for the case of ideal gas can be 
interpreted as an effective dimensional reduction of the 
space at the high temperature regime which is a com¬ 
mon feature of quantum gravity proposals [Dl] and also 
phenomenological approaches to the minimal length con¬ 
jecture [47]. 

Considering the limit 0 = 1 is useful at least 

for two reasons. First, we see that the modified thermo¬ 
dynamical relations recover the standard results at the 
sufficiently low temperature regime when the minimal 
length effect are negligible (correspondence principle). 
Second, one can estimate the magnitude of the order of 
the quantum gravity corrections to the thermodynami¬ 
cal quantities which may potentially lead to observable 
effects. In this limit, the series in the relation (l43l) con¬ 
verges and one can use it to obtain all the minimal length 
(quantum gravity) corrections to the internal energy and 
specific heat. Substituting the relation (H51) into Cfil) and 
then using the definition U = T 2 ( dl ! ^ r Z ) N , it is easy to 
show that the internal energy in the limit of 0 ^ 1 will 
be 



n(2n + 1)!! 
( 20 2 )” ’ 


(47) 


where Uq = 3NT is the internal energy for the standard 
nondeformed N independent three-dimensional harmonic 
oscillators. In this limit, the specific heat can be also 
obtained by substituting the above relation into the def¬ 


inition C v = (§ y) v which gives 


Cy 2 , n(n + l)(2n + 1)!! 

Coy “ 3 ^ j (20T 


(48) 


where Cq v = rOy = 3-/V i s the specific heat of the 
corresponding nondeformed case. 

As is clear from the relations (l47l) and (l48l) . the first 
minimal length (quantum gravity) corrections to the in¬ 
ternal energy and specific heat of the harmonic oscillator 
are of the order of 0 -2 . To estimate the magnitude of 
the order of these corrections, consider the vibrational os¬ 
cillations of a carbon monoxide molecule which may be 
modeled with the oscillators of mass m ~ 10~ 26 kg and 
frequency cr ss 10 15 Hz. With this numerical value for 
the mass and also considering k = 0( 1 )T P1 ~ 10 19 GeV 
and T ~ lTeV, we get 0 ~ 10 17 . Therefore, the first 
minimal length corrections to the internal energy and 
specific heat of the harmonic oscillator in Snyder space 
are of the order of 0~ 2 ~ 10 -34 which is too small to be 
experimentally detected for the accessible energy scales. 


V. SUMMARY 

Appearance of the Hamiltonian systems with nontriv¬ 
ial structure in the context of phenomenological quantum 
gravity candidates, such as the doubly special relativity 
theories with deformed noncommutative phase spaces, 
naively suggests the revision of the statistical mechan¬ 
ics formalism. In the first step, we have formulated the 
statistical mechanics of a general Hamiltonian system 
in a covariant (chart-independent) manner by means of 
the symplectic geometry. The results show that the two 
properties of the phase space, as a symplectic manifold, 
play distinguished roles in statistical consideration of a 
system: The topology of the phase space as a global prop¬ 
erty and the local form of the symplectic structure that 
determines the geometry of the phase space. For a topo¬ 
logically trivial phase space with the canonical symplectic 
structure, the standard statistical mechanics emerges as 
it is expected. The subtleties arise when the phase space 
has nontrivial topology or geometry such as the phase 
spaces with curved momentum space and noncommuta¬ 
tive structure. The topology of the phase space turns 
out to be related to the total number of microstates, for 
instance, the number of accessible microstates will be 
finite for a system with compact phase space. The sym¬ 
plectic structure also affects the probability distribution 
of the microstates. The canonical form of the symplectic 
structure in terms of the positions and momenta of the 
particles plays an important role in the standard statisti¬ 
cal mechanics. It induces a uniform probability distribu¬ 
tion over the microstates by defining a uniform measure 
on the corresponding phase space. However, the non¬ 
commutative phase space always provides a noncanon- 
ical symplectic structure in terms of the positions and 
momenta of the particles and the probability distribu- 
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tion then will be nonuniform. This is a general result 
for the phase spaces with noncanonical (noncommuta- 
tive) structure which can be immediately realized from 
the covariant formulation of the phase space. We imple¬ 
mented the covariant formalism in order to study the 
statistical mechanics in Snyder noncommutative space 
as an explicit example of a phase space with nontriv¬ 
ial topology and geometry. We obtained the associated 
partition function from which all the thermodynamical 
properties of the system can be obtained. As a particular 
example, we obtained the partition function for the three- 
dimensional harmonic oscillator and we have shown that 
our result is in good agreement with that which arises 
from the full quantum consideration at the high tem¬ 
perature regime. While there is no analytical expression 
for the full quantum partition function of the harmonic 
oscillator in Snyder space, the semiclassical approxima¬ 
tion provides an analytical partition function which is 
applicable at the high temperature regime. Using the 
obtained partition function, we have studied the ther¬ 
modynamical properties of the system of harmonic oscil¬ 
lators in Snyder space and our analysis shows that the 
number of microstates will be drastically reduced in this 
setup due to the existence of a minimal length. This re¬ 
sult justifies our general claim for the deformed spaces 
which take into account a minimal length scale. Apart 
from the details of the models which deal with a univer¬ 
sal minimal length such as the noncommutative spaces, 
doubly special relativity theories, the generalized uncer¬ 
tainty principle, and polymer quantization, the covariant 
formalism reveals the main role of the minimal length in 
statistical mechanics: The microstates with higher en¬ 
ergy or momenta are less probable when there is a mini¬ 
mal length scale for the system. Also, we calculated the 
quantum gravity corrections to the internal energy and 
specific heat of a system of harmonic oscillators. We esti¬ 
mated the order of magnitude of the first minimal length 
corrections which are of the order of 10 _34 for both of 
the internal energy and specific heat. The equilibrium 
as an essential criterion for the statistical system is also 
considered in the Appendix by means of the dynamically 
invariant Liouville volume which enters in the definition 
of the Gibbs entropy and partition function of the statis¬ 
tical system. 


APPENDIX: EQUILIBRIUM 
A. Statistical independence 

Relation (l28l) can be seen as a particular case in which 
the system is divided into statically independent sub¬ 
systems and subsystems were taken to be the particles 
themselves. It is clear that the statistical independence 
between particles will be spoiled for a system in which 
the particles interact with each other. However, when 
the system is in equilibrium, one can still divide it into 
M macroscopic subsystems (1 -C M <C N) which are 


statistically independent. In this respect, the total phase 
space T and the associated symplectic structure w will be 

M 

r = ri x ... x r M , u = (a-i) 

A—l 

The corresponding Liouville volume is then given by 

uj 3N = (wi) 3 " 1 A ... A( w „) 3 "“, (A-2) 

where n A is the number of particles of the subsystems and 
is assumed to be sufficiently large (n A 1) to guaran¬ 
tee that the subsystems are macroscopic. The statistical 
independence of the subsystems is then defined as 

pco 3N = p 1 (wi) 3ni A ... A p M ( oj M ) 3nM ■ (A-3) 

This criterion is usually employed as an intrinsic charac¬ 
terization of an equilibrium probability distribution 0] . 
Indeed, the system is in equilibrium if components of the 
system are in relative equilibrium with respect to each 
other. Based on this internal definition of the equilib¬ 
rium, the Gibbs state can be understood as an equilib¬ 
rium state as follows. Consider the total Hamiltonian 
function of the system as 

M 

H=J2h A , (A-4) 

A=1 

where {Ha} are the Hamiltonian functions of the M 
macroscopic subsystems and we also disregard the in¬ 
teraction between subsystems based on the fact that the 
relative number of particles which take part in the in¬ 
teractions is negligible compared to n A . So, using the 
relations ESI) . (IA-1D . and (IA-41) . the Gibbs state turns 
out to be an equilibrium state based on (IA-31) . Addition¬ 
ally, it is straightforward to show that 

[x£,xf]=0, (A-5) 

which is the criterion for the kinematical and dynami¬ 
cal separability of the subsystems. The relation (IA-5D 
then can be viewed as the equilibrium criterion when the 
system is in a Gibbs state. 

Note that the decomposition (IA-41) holds only over not- 
too-long intervals of time since the effects of interaction of 
the subsystems will eventually be dominated even if such 
interactions are too weak such as the universal gravita¬ 
tional effects 0. 

B. Liouville equation 

Another more common definition for the equilibrium 
statistical state of the system is 



in which p e denotes the density that corresponds to the 
equilibrium state. The fact that the Gibbs state is an 
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equilibrium state then will be justified through the Liou- 
ville equation 

d ( ^ 3 ») = (| + £,„)(^ 3 '*)=0, (A-7) 

in which the Liouville theorem m is used. The Liou¬ 
ville equation should be satisfied by any statistical state 
as well as Gibbs one and it is, indeed, the necessary re¬ 
quirement for the consistent probabilistic interpretation 
of p. Using the fact that £ x p = x. H p = {p,H}, the 
Liouville equation (IA-71) can be rewritten in the more 
well-known form 

| + {p,tf} = 0. (A-8) 

From the relation (E6l) . it is clear that {p c ,H} = 0, 
and the Gibbs state then satisfies the relation dp c /dt = 
0, which shows that it is indeed an equilibrium state 
through the definition (IA-61) . 


If we use the Liouville equation (IA-71) and consider 
the fact that the Liouville volume oj 3N does not change 
explicitly with time, it turns out that 


dS_ 

dt 



!+£ x „) J^p\np uj 3N 
P C Xh co 3N = 0, 


(A-9) 


which shows that the fine-grained entropy (l?5l) does not 
change with time as expected for closed systems. The 
above relation and the Liouville equation (IA-71) guarantee 
that the number of microstates does not change through 
the dynamical evolution of the system, which makes the 
statistical formulation a consistent setup to give the ther¬ 
modynamical properties in any ensemble. 
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